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Abstract 

The difference in production rate between and W~ at hadron col- 
liders is very sensitive to the the difference between up- and down-quark 
distributions in the proton. This sensitivity allows for a variety of use- 
ful measurements. We consider the difference Q^) — Us{x,Q'^) in the 
sea distributions and the difference /S.u{x, Q"^) — Ad{x, Q"^) in the polarized 
parton distribution functions. In both cases we construct an asymmetry to 
reduce systematic uncertainties. Although we discuss measurements at the 
Tevatron and future hadron colliders, we find that the Brookhaven Rela- 
tivistic Heavy Ion Collider (RHIC) is the most appropriate hadron collider 
for these measurements. 
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1 Introduction and Results 



It has been appreciated for some time that useful information on the parton distri- 
butions of nucleons can be extracted from measurements of the charge asymmetry in 
hadronic W production |1|, ^, ^ The asymmetry provides us with information on 
the difference between the up- and down-quark structure of the proton, i.e. on the 
quantities 6qy = — d^ and 6qs = dg—Us- Here, as usual, the subscripts v and s refer 
to the valence and sea components of the proton structure functions. The quantity 6qy 
is studied in deep inelastic scattering experiments on proton and neutron targets and 
is critical to the determination of the hadronic width of the W from the ratio of the 
W and Z cross sections in hadron collisions. While it is often assumed that Sqs = 
on the basis of isospin symmetry, this does not have to be the case. As ^ d^, 
evolution will inevitably result in Us ^ dg. Also, as there are more u- than (i-quarks in 
the proton, one might imagine that the further generation of uu pairs is suppressed by 
Fermi statistics]^]. As we discuss further on, experimental evidence for the violation 
of the Gottfried sum rule has been interpreted as being the result of the non-vanishing 
of 

Analogously for W production with polarized beams, a double asymmetry in the 
initial proton helicities and in the W charge can provide us with useful information 
on the up- and down-quark structure in the polarized proton. One of the theoretical 
assumptions generally used in the interpretation on the EMC result on proton spin 
is based on a combination of low energy experiment and isospin symmetry, i.e. the 
difference between the polarized u- and d quark distributions. This quantity can be 
directly measured at higher energy in W production at a polarized proton-proton 
collider. 
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The main point of this paper is that the W charge asymmetry in hadron colhsions 
measures both Sq^ and 6qs. We show, moreover, that the latter quantity can only 
be effectively probed in proton-proton colliders. W production at RHIC thus has 
unique aspects. This includes the possibility of using polarized beams, which we discuss 
in some detail. 



2 Gottfried sum rule in unpolarized production 

The recent NMC0] result on the Gottfried p|] sum rule 



(extrapolated to x = 0) of /g(0) = 0.240 ± 0.016 is significantly different from the 
result of 1/3 expected if the usual assumptions are made, i.e. the proton consists 
of valence uud and u{y) = d{y). This "violation" of the Gottfried sum rule is quite 
easy to understand — there are more w-quarks in the proton than rf-quarks, so the 
Pauli exclusion principle suppresses the production of uu pairs in the sea relative to 
dd pairs]^ — but an independent confirmation of this result is desirable. 

The large number of W's (we estimate 3000 for = 200 GeV and 300,000 
for ^/s = 500 GeV) to be produced at the Brookhaven Relativistic Heavy Ion Col- 
lider (RHIC) will provide the necessary confirmation as has been suggested in the 
literature^. Instead of considering the separate production of W'^ and W~, consider 
the charge asymmetry v4iy[|[]: 




(1) 



Dw 



Nw 



Aw 



da{W+) _ da{W-) 

dyw dyw 
da{W+) dajW-) 

dyw dyw 
Nw 
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It is straightforward to show that (to simphfy the notation we suppress the dependence 
of the parton distribution functions on in the following): 

Nw OC COi^^ dc{Us{x2) [Uy{xi) - d^{xi)\ + Uy{x2) [ds{xi) - Us{Xi)]} 
+ SIV? 9c{Uy{Xi)Ss{x2) - Cs{xi)d.^{x2)] 

+ (xi ^ X2) (3) 

and 

Dw OC cos'^ 9c{Us{x2) [uv{xi) + d^{xi)] + Uv{x2) [ds{xi) + Us{xi)] + 2Ss{xi)Cs{x2)} 

+ SW? ec{u{Xi)Ss{x2) + Cs{Xi)d{x2)} 

+ ^ X2), (4) 

where 6'c is the Cabibo angle. This gives the exact expression we use to calculate Aw- 
In calculating these cross sections we include the usual i^'-factor for W production, 
K = 1 + 87ra,(M2,)/9 (see, e.g., Ref. fTp). 

To further illustrate our analysis, consider only the terms proportional to cos^ Oc- 
This is a relatively good approximation in any case as sin^ Oc ~ 0.05. In this limit, the 
numerator factor reduces to 

Nw ~ Us{x2)5q.u{xi) + u^{x2)5qs{xi) + (xi ^ X2) (5) 



^ Mw 



where dq^i^x) = Uy{x) — dy{x) and 5qs{x) = ds{x) — Us{x). W production is central and 
relatively flat in rapidity, so we examine Aw at y = 0. Here, xi = X2 = xq 
and 

^ I _ Us{xo)5qy{xo) + Uy{xo)Sqs{xo) 
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A similar asymmetry can be formed at pp colliders. In this case, the numerator 
factor is 

Nw oc cos^ 9c{us{x2)dq^,{xi) + ^Sqv{xi) {u^{x2) + dy{x2)) + Uy{xi)Sqs{x2)} 
+ siri^ 9c{uy{xi)ss{x2) + Cs{xi)dy{x2)} 

- {xi ^ X2). (7) 

The expression separates into terms proportional to Sq^, and 6qs and another small term 
proportional to sin^ 9c, but it is not as simple as the expression (5) for the asymmetry 
at pp colliders. More importantly, the asymmetry vanishes for y = at pp colliders. 
In the central region of rapidity, where most of the W^s are produced, the asymmetry 
is small, severely limiting the usefulness of pp colliders in this measurement. 

We have demonstrated the usefulness of pp coUiders in the measurement of the 
difference in the distributions Us and ds- RHIC has been designed with a variable 
center-of-mass energy, ranging from 50-500 GeV. This gives an xq range 0.16 < xq < 
1 (though with low statistics at threshhold). SSC and LHC are sensitive to Xq — 
0.002 and 0.005 respectively (although moving away from y = allows for a range 
of sensitivity about these values). At such small values of x, sea quarks dominate W 
production. Referring to Eqn. (6), it is clear that the term Us{xo)5qy{xo) will easily 
dominate the asymmetry, leading to very little sensitivity to Sqs{xo) at SSC and LHC. 
This argument does not hold at RHIC since the larger x region probed leads to a 
larger valence contribution, which in turn leads to a rather large contribution from the 
difference in the light sea quark distributions. Thus RHIC is the only pp collider where 
this measurement can be performed. 

Our conclusion thus far is that RHIC is uniquely suited for this measurement, 
because a) Aw is large at small rapidity where the statistics are best, unlike all pp 
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colliders, b) RHIC probes the correct x-range, unlike pp supercolliders and c) RHIC 
has variable center-of-mass energy, which allows for measurement at different x without 
measuring at different rapidity. 

3 Results on the Gottfried Sum Rule 

Consider the asymmetry Aw at the Tevatron. In Fig. 1 we show the cross section for 
and W~ production, for various parton distribution functions. Next, we show the 
numerator factor N\y and the dominant contribution from 6qs (Fig. 2). The contribu- 
tion from a non-5'f/(2) symmetric sea is very small compared to the contribution from 
the difference in valence distributions, As we need to measure this asymmetry at 
large rapidity (~ 1.5), we see that information on dg — Ug will be difficult to obtain 
from the Tevatron. 



Some data exist on the leptonic asymmetry at the Tevatron. Here, 



N, 



dye dye 
da{i-) da{£+) 



A, 



dye dye 
Ne_ 
De 



is the charge asymmetry in rapidity distribution for single lepton production. We show 
in Fig. 3 our result for Ag, the asymmetry for electrons and compare it to the data 



from CDF||TO|. In Fig. 4 we compare the contribution to N^. from 5qs with the total 
numerator. 

Next, consider the asymmetry Aw at hadron supercolliders. Our results for the 
SSC are also representative of the results that can be obtained at the LHC In Fig. 5 
we show the numerator factor Nw and the dominant contribution from 5qs. It is clear 
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that the contribution from a non-SU{2) symmetric sea is again very small compared 
to the contribution from the difference in valence distributions, Although the 

statistics will be very good, the relative smallness of 6qs will make this measurement 
difficult. 

Finally, we consider the asymmetry Aw at RHIC, a high luminosity [C = 2 x 
10^^ cm~^sec~^ = 6000 pb^^/yr) collider capable of producing proton-proton collisions 
for center-of-mass energies between 50 and 500 GeV. RHIC will also collide protons 
on heavy nuclei in order to connect to existing data and heavy ions on heavy ions. 
We assume a nominal running time of two months at full luminosity. In order to be 
somewhat conservative, we estimate event numbers based on 300 pb~^ integrated lu- 
minosity. We assume a generic collider-type detector, and we require the photons and 
electrons observed to lie in the rapidity range \y\ < 2. This simulates the acceptance 



of the proposed STAR detector at RHIC, level 2 for photons and electrons [|T2[. We do 
not consider the possibility of the detection of muons at RHIC. We show the numer- 
ator factor Nw, and the contribution to from 6qs in Figs. 6a and 6b respectively 
for ^/s = 200 GeV and 500 GeV. At RHIC, the contributions from 6qs are in no way 
suppressed; they can in fact dominate Nw- The parton distributions with a symmetric 
sea {e.g., HMRSB, EHLQ2 and MRSDO) generally give smaller asymmetries (espe- 
cially at small W rapidity, yw) than the distributions with non-symmetric sea {e.g., 
MRSDM, CTEQIL and CTEQIM). Fig. 7 gives the asymmetry Aw for = 500 GeV 
for various choices of the parton distribution functions. Here at last we can observe 
in the asymmetry the effects of the non-symmetric sea, although for this particular 
scenario, the 6qs contribution using the CTEQIL distributions accidentally small. The 



explanation for this phenomenon is given in Ref. - the sign oi d — u changes with 



X, and W production at ^/s = 500 GeV is sensitive to values of x near the crossover 
point for the CTEQIL distributions. 

Having determined that this measurement is in principle possible at RHIC, we study 
events with a single high-p^, lepton and missing p^. In Figs. 8a and 8b we show the 
rapidity distributions for production of e"*" and e~ at center of mass energies of 200 and 
500 GeV, respectively. Even in the worst case (e~ production at 200 GeV), we expect 
about 15 events in the rapidity range \ye\ < 0.5, so the statistics should be sufficient. 
In Figs. 9 we compare the contributions to N^. Again, the parton distributions with a 
symmetric sea give generally smaller asymmetries that the parton distributions with 
non-symmetric sea, as can be seen from the Sqs contributions in Fig. 9, although the 
accidentally small Sqs contribution from CTEQIL distributions at ^/s — 500 GeV 
is apparent. Finally, in Fig. 10 wc show the observable asymmetry Ae for center- 
of-mass energies of 200 and 500 GeV. Here, especially at small electron rapidity ye, 
the asymmetries from the symmetric sea parton distributions cluster together and 
are reasonably separated from the asymmetries from the non-symmetric sea parton 
distributions. 

4 Au — Ad in polarized production 

The interpretation of the EMC result on proton spin relies upon two low energy quanti- 
ties, namely (Am— Ao?) = (from isospin invariance) and {Au+Ad—2As) = 3 F — D 
(from SU (3) symmetry) where qa is the axial- vector coupling in neutron /3-decay and F 
and D are the invariant amplitudes for the axial- vector current in hyperon semileptonic 
decays. Here, in the usual notation 

Aq{i^'^)sa = {p,s\ q'^albq | P, s) ^2, (9) 
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where | p, s) is a proton state with spin vector Sa- The quantity Aq is related to the 
polarized parton distribution function for quark q 

Ag(g') = ['Aq{x,Q^)dx (10) 
Jo 

where Aq{x, Q^) is the difference in the distributions of a quark q in a longitudinally 
polarized proton with the same and opposite helicity as the proton. As long as the 
assumptions of isospin invariance and SU (3) symmetry are valid, these particular (non- 
singlet) combinations of the Aq''s will not run with Q'^; they reduce to the difference 
and sum of the first moments of the polarized up- and down-valence distributions. 
Thus, the two quantities above can be used at higher energy even thought they are 
extracted from low energy data. It would, nonetheless, be comforting to extract one 
(or both) of these quantities at a higher energy as a consistency check. 

Analogous to the earlier discussion of unpolarized W production we can derive the 
following for polarized W production: 

(X cos^ 9cAu{xi)Aq{x2) + sin^ 9cAu{xi)As{x2) + {xi ^ X2) 

(11) 

oc cos^ 9cAd{xi)Aq(x2) + sin^ 9c Aq{xi)As{x2) + {xi ^ X2)- 

(12) 

We assume that the charm content of the proton is zero, which is a very good approx- 
imation below Tevatron energies, and Au = Ad = Aug = Adg = Aq, and wc again 
suppress the dependence of the parton distriubtion functions on Q^. In the usual 
notation, a'^'^ (c""*"") is the cross section for two protons with the same (opposite) he- 
licit ies and the polarized parton distribution functions as described above. As before, 
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da++{W+) _ da+-{W+) 
dyw dyw 

da++(W-) _ da+-(W-) 
dyw dyw 



we include the standard X-factor in our calculation of the cross section. Just as the 
difference of cross sections for unpolarized W production is the most interesting and 
useful quantity, we find that the difference 

da++{W+) _ da+-{W+) 
dyw dyw 



da++{W-) da+-{W-) 1 
dyw dyw J 



{ cos^6lc (A'u(xi) - Ad{xi)) Aq{x2) + sin^ 6lc Aii^(xi) As{x2) + {xi ^ xaX^S) 

is particularly useful in isolating the differences we would like to measure. 
Consider the quantity All- 



All = 



da++{W+) _ da+-{W+) 



da++{W-) _ da+-{W- 
dJj^ dyw 



da++{]V+) , da+-(]V+) ^ da+-{W-) , da+-{]V-) 



(14) 



dyw ^ ^ dy^ ^ ^ 

This asymmetry is again dominated by the cos^ 9c terms, not only because sin^ 9c is 
small but also because Au{x) and Ad{x) have opposite signs, so the interference is 
constructive. In the limit that we only consider the cos^ 9c terms, this asymmetry 
reduces to 



A 



LL 



[Au{xi) — Ad{xi)]Aq{x2) + {xi ^ X2) 



(15) 



[u{xi) + d{xi)]q{x2) + {xi ^ X2) ]' 
Finally, because the cross section for central Ws is large, we can examine this asym- 
metry at y = 0, which gives: 

[Au{xq) - Ad(a;o)]Ag(xo) 



A 



LL \y=Q 



[u{xq) + d{xQ)\q{xQ) 



(16) 



where a;o = ^ as before. Inclusion of the sin 9c terms will slightly modify this 
result, but it should be possible to extract the combination Au{x) — Ad{x) from the 
data due to the large statistics. Furthermore, by running RHIC at several center-of- 
mass energies and moving away from y = 0, it will be possible to map the desired 
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polarized parton distributions for a range of x values. We give the expressions above 
under the approximation sin^ 6c = merely to simplify the discussion. In the figures 
we use the exact expressions. 

5 Results on Au — Ad 

A program of polarized proton-proton collisions, at full energy and luminosity, is being 



discussed |1TT|]. We show in Fig. 11, the W rapidity distribution, da/dyw, for and 
W~ production at two representative RHIC energies, 500 GeV and 200 GeV. Given 
the large integrated luminosity from a two month run at RHIC (~ 300 pb^^), there 
will be a rather large number of W^s produced at RHIC. Next, we present our results 
for All as a function of the W rapidity yw for RHIC at y/s = 200 GeV (Fig. 12a) and 
at ^/s = 500 GeV (Fig. 12b). In Fig. 12, we use two extreme cases for the polarized 
parton distribution functions, one in which the polarized gluon is taken to be large (and 
consequently the polarized strange quark distribution is small - set ls=0 of Bourrely, 
Guillet and Chiappetta (5(7(7) |]TB[) and one in which the polarized gluon is small 
(and the polarized strange quark distribution large - set 1^=1 of BGC). We note that 
the large difference in the asymmetries calculated using the different polarized parton 
distribution functions is due primarily to the large uncertainty in the polarized sea 
distribution. This uncertainty can be traced back to the different explanations of the 
EMC effect on proton spin. This uncertainty will be reduced using other measurements 
at RHIC, e.g., jetfl^, |1^, 0, direct photon [p!7|, |18|, |19|, heavy quark [pO[l or charmonium 



production ||2T|, |2^, designed to measure the polarized gluon contribution to 



proton spin. Once this error is reduced, it will be possible to extract the desired 
quantity, Au — Ad. Also, the polarized parton distribution set with the 'smaller' 
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polarized sea quark gives a larger asymmetry at large rapidity for y/s = 200 GeV. This 
is a consequence of the different large-x behavior in the polarized parton distributions; 
the event rate at large rapidity is sufficiently small that for practical purposes, this 
region of large rapidity will contribute very little to the observed asymmetry. 

Of course, the experiments at RHIC will not directly observe W^s, but rather will 
reconstruct them from the sample of single isolated lepton and missing p^, so one should 
also study the possibilities including the decay. To simulate the aceptance of the STAR 
detector, we impose a cut of |?/| < 2 on electrons; STAR will have no acceptance for 



muons 



. In addition, we place a cut on the transverse mass of the reconstructed W 
of 50 GeV, which will remove much of the background without significantly decreasing 
the signal. In Fig. 13 we show the rapidity distribution, da/di/e, for single electrons 
and positrons for two representative center of mass energies, 500 and 200 GeV. Even 
including the W decay there will still be a large number of events (about 130 for 
= 200 GeV and 13000 for ^/s = 500 GeV). We give, in Fig. 14 the asymmetry 
Ali^ including electronic decay of the W, with = 200 GeV (Fig. 14a) and with 



^/s = 500 GeV (Fig. 14b). The curves follow the convention of Fig. 12. Here, 

(17) 



A 



dff++(e+) rfa+-(e+) 
dye dye 



rf(T++(e-) _ da+-{e-) 
dye dye 



1 da++{e+) da+'{e+) da++{e-) da^-^e- 

'(Efe We We We 

Again, the polarized parton distributions that give use a large polarized gluon to 
explain the EMC effect on the spin of the proton give give a small asymmetry due 
to the presence of a Ag factor in the asymmetry, and at \/s = 200 GeV the different 
large-x behavior of the polarized parton distributions is apparent in the crossings of 
the asymmetries. Once the uncertainties in the polarized gluon and sea distributions 
are reduced, it will be possible to extract Am — Ad. 
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Figure Captions 

Figure 1 Rapidity distribution for the production of W" and at the Tevatron. The 
different curves correspond to different parton distribution functions: dotted 
- HMRSB[|1], double dotted - EHLQ set 2|2§, sohd - MRSDO§, dashed - 



MRSDMfl, dot-dashed - CTEQ1L[|3] and double dot-dashed - CTEQ1M[|7 



Figure 2 Numerator factor Nw (as defined in the text) versus yw at the Tevatron. The 
different curves follow the convention of Figure 1. Both the total N\y and the 
contribution from 6qs are shown. 

Figure 3 Electron asymmetry at the Tevatron. The labeling of the different curves 
follows the convention of Figure 1. Data are from Ref. |[10|| . 



Figure 4 Numerator factor Ne (as defined in the text) versus t/e at the Tevatron. The 
different curves follow the convention of Figure 1. Both the total numerator 
factor and the contribution from 6qs are shown. 

Figure 5 As Fig. 2, at SSC. 

Figure 6 As Figs. 2, 5, for a) RHIC at 200 GeV and b) RHIC at 500 GeV. 

Figure 7 The charge asymmetry Aw (as defined in the text) versus yw at RHIC for i/i = 
500 GeV. The different curves follow the convention of Figure 1. 

Figure 8 The rapidity distribution for single e"*" and production at RHIC. The dif- 
ferent curves follow the convention of Figure 1. Figure 8a shows the rapidity 
distribution for y/s = 200 GeV and Figure 8b show the rapidity distribution for 
= 500 GeV. 
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Figure 9 As Fig. 4, for a) RHIC at 200 GeV and b) RHIC at 500 GeV. 

Figure 10 The charge asymmetry Ag in single lepton production (as defined in the text) 
versus i/e at RHIC. The different curves follow the convention of Figure 1. Fig- 
ure 10a shows Ae at i/i = 200 GeV and Figure 10b shows Ag at i/i = 500 GeV. 

Figure 11 The rapidity distribution for production of W~ and at RHIC, using EHLQ, 



set 2^6[ parton distribution functions. The solid curve corresponds to pro- 
duction at y/s = 500 GeV, the dashed curve corresponds to W~ production at 
■\/s = 500 GeV, the dotted curve corresponds to iy+ production at ^/s = 200 GeV 
and the dot-dashed curve corresponds to W~ production at y/s = 200 GeV. 

Figure 12 The asymmetry All (as defined in the text) versus yw at RHIC. Figure 12a 
shows All for = 200 GeV and Figure 12b shows All for ^/s = 500 GeV. 
We use EHLQ, set 2|^ for the unpolarized parton distributions, while in both 
cases, the solid curve uses the polarized parton distribution functions of BGC, 
set Is = and the dashed curve used the polarized parton distribution functions 
of BGC, set Is = l- 

Figure 13 The rapidity distribution for production of e~ and e"*" at RHIC, using EHLQ, set 



2[^ parton distribution functions. The solid curve corresponds to production 
at = 500 GeV, the dashed curve corresponds to e~ production at ^/s = 
500 GeV, the dotted curve corresponds to e+ production at ^/s = 200 GeV and 
the dot-dashed curve corresponds to e~ production at y/s = 200 GeV. 

Figure 14 The asymmetry Aj^j^ (as defined in the text) versus i/e at RHIC. Figure 14a shows 
AIl for ^/s = 200 GeV and Figure 14b shows AIl for ^/s = 500 GeV. We use 



EHLQ, set 2|26| for the unpolarized parton distributions, while in both cases, the 
sohd curve uses the polarized parton distribution functions of BGC, set Is = 
and the dashed curve used the polarized parton distribution functions of BGC, 
set Is = 1. 
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